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In this paper we review some results of our works on the magnetization processes in: i) Isolated 
nanomagnets, both in the one-spin approximation and as many-spin systems. Here, we focus on 
the intrinsic properties, e.g., those induced by finite-size, boundary and surface effects. We also 
investigate the crossover between the two regimes, ii) Assemblies of nanomagnets, also in the two 
situations. We focus on their behavior mainly due to dipole-dipole interactions. Then, we will 
comment on the interplay between these intrinsic and collective effects. 
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I. INTRODUCTION 

Magnetic information storage is one of the major ap- 
pUcations of magnetic nanoparticles, though magnetic 
nanoscalc structures in general have been identified as 
offering a great potential in appUcations ranging from 
nano spin electronics, innovative complex materials like 
magnetically controlled shape memory alloys to medi- 
cal applications in medical diagnosis and therapy. The 
objective of achieving high-density requires rather small 
sizes, high anisotropics, and imposes severe constraints 
on the microstructurc and the magnetic properties to 
ensure a good stability against thermal demagnetization 
over several years. On the other hand, the same material 
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must be magnetized by the writing head in a time inter- 
val as short as a nanosecond. However, smaller sizes in- 
duce stronger surface effects and higher densities induce 
stronger inter-particle interactions. Unfortunately, sur- 
face effects and inter-particle interactions have the draw- 
back to negatively affect the signal-to-noise ratio. Opti- 
mising the underlying material in view of technological 
applications requires taking account of and understand- 
ing the effects of surface and interactions, both on the 
static and dynamic level. 

At present there are mainly two types of nanoparti- 
clc samples: i) Assemblies of nanoparticles embedded 
in a non-magnetic matrix with volume distribution and 
randomly oriented easy axes, with negligible to strong 
dipole-dipole inter-particle interactions (DDI) Q . ii) Iso- 
lated single particles measured by the technique of /i- 
SQUID Q. From a theoretical point of view, none of 
these sample prototypes is so far fully satisfactory: While 
in the case i) one can measure by commercial SQUID 
most of the magnetic properties, e.g., field-cooled (FC), 
zcro-field-cooled (ZFC) magnetization and ac suscepti- 
bility, in a wide temperature range (4.2-300K) and in 
any applied magnetic field available (0-30T), this case 
presents tremendous difficulties for modeling because of 
the volume and anisotropy easy axes distributions and 
DDI. Indeed, these interactions, because of their long- 
range character, together with the two distributions, re- 
quire prohibitively long CPU times to deal with realistic 
assemblies of, say 10"' particles of 300 — 1000 spins each. 
Besides, the physical quantities, such as the magnetiza- 
tion of the assembly, are only known on average. In the 
case ii), one gets rid of these difhculties but then the 
experimental results are limited to low temperature and 
applied field, and more importantly the magnetization 
itself is not accessible, and for very small particles the 
signal is rather "noisy" , which hinders any measurement 
of hysteresis loops, and one makes do with the angular 
dependence of the switching field, the so-called Stoner- 
Wohlfarth astroid. 

In many theoretical approaches the magnetic state of 
a small magnetic particle is represented by a single mag- 
netic moment regardless of its microscopic origin, this is 
the one-spin approximation^ or one-spin problem (OSP), 
which can only be justified for particles that arc not 
too large to be single-domain, and not too small to be 
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surface-effects free. Within this approximation, it is well 
known by now that the magnetization of a nanoparticle 
can overcome the anisotropy-energy barrier and thus re- 
verse its direction, at least in two ways: either under ap- 
plied magnetic field which suppresses the energy barrier, 
or through thermally activated statistical fluctuations. 
Switching under applied field, at zero temperaturejis de- 
scribed by the Stoncr-Wohlfarth (SW) model [1, In 
fact, it has recently been shown experimentally and 

theoretically, that efficient magnetization switching can 
also be triggered by transverse field pulses of a duration 
that is half the precession period. The SW model ac- 
counts for the hysteretic reversal of a magnetic moment 
over the potential barrier under the influence of the mag- 
netic field applied in an arbitrary direction and at zero 
temperature. In fact, at zero temperature the magneti- 
zation switching is only possible when the field reaches 
a critical value that suppresses the energy barrier. At 
finite temperature and short-time scales, crossing of the 
energy barrier is triggered by statistical thermal fluctu- 
ations and is described by the Neel-Brown (NB) model 
0, H, H, and its extensions, which took several years 
to be fully accomplished [see the review article and 
references therein]. Both the SW and NB models have 
been confirmed, to some extent, by experiments on indi- 
vidual cobalt particles [2| through measurements of the 
switching field distribution, the non-switching probabil- 
ity, and the time evolution of the direction of the parti- 
cle's magnetization. At finite temperature, but at quasi- 
equilibrium, the magnetization switching occurs accord- 
ing to two distinct regimes At very low tempera- 
ture, switching operates by coherent rotation of all spins, 
as in the SW model, whereas at higher temperatures, the 
magnetization switches by changing its magnitude. The 
latter effect results in a shrinking of the SW astroid as 
qualitatively described by the modified Landau theory 
and confirmed by experiments 

However, it is clear that the change of the magnetiza- 
tion magnitude during its switching cannot be explained 
in the framework of the OSP approximation. Indeed, 
it can only be understood as the result of a successive 
switchin g o f individual (or clusters of) spins inside the 
particle |l3l . [13] , which should then be necessarily dealt 
with as a many-spin system. In more general terms, the 
picture of a single-domain magnetic particle with all spins 
pointing into the same direction, leading to coherent re- 
laxation processes becomes unsatisfactory as soon as one 
comes to deal with small particles with a surface that 
makes up to 50% of the total volume. This is to say, the 
surface entails strong effects that cannot be neglected. 
In fact, deviations from the OSP approximation, and 
thereby from both the SW and NB models, have been 
observed in metallic particles [isl. [T6| . ferrite particles 
(itI . [ist , and maghemite particles [l9l ]. These deviations 
have materialized in terms of the absence of magneti- 
zation saturation at high fields, shifted hysteresis loops 
after cooling in field, and enhancement at low temper- 
ature of the magnetization as a function of the applied 



field. 

In most cases, surface effects happen to be strong 
enough as to compensate for the work needed against 
the exchange energy that favors full alignment, and it 
is conceivable then that the magnetization vector points 
along the easy axis in the core of the particle, and then 
gradually turns into a non-coUinear direction when it ap- 
proaches the surface. Moreover, in addition to exchange 
interactions, which are the strongest interactions between 
atomic moments in a magnetic system, there are also 
the purely magnetic dipole interactions between the mag- 
netic moments of the atoms and the interactions between 
the magnetic moments and the electric field of the crys- 
tal lattice (spin-orbit interactions). The last two types of 
interactions, being relativistic in origin, are much weaker 
than the exchange energy on a short-distance scale. How- 
ever, on a long-distance scale they are non-negligible be- 
cause they are of long range. These interactions have two 
important features: they induce inhomogencities in the 
spatial distribution of the magnetization, and introduce 
a preferred direction in the system Roughly, DDI 

lead to two energy terms corresponding to the volume 
and surface charges. In a small magnetic system, such as 
a nanoparticle, only the second contribution is significant 
and it accounts for the shape anisotropy, and the former 
becomes negligible as one integrates over a small volume 
[21I . [2^ . Therefore, to study a small magnetic system 
where inhomogeneities are necessarily present, one has 
to take into account all different contributions to the en- 
ergy: exchange and DDI, core and surface anisotropics, 
and Zeeman energy. 

So, one of our goals is to understand surface effects on 
the thermodynamic and spatial behavior of the magneti- 
zation in small systems. This is also of crucial importance 
to the study of their dynamics. However, this requires 
a microscopic approach that accounts for the local en- 
vironment inside the particle and the above mentioned 
contributions. Unfortunately, this leads to a rather diffi- 
cult, if not insuperable, task owing to the large number 
of degrees of freedom which hinders any attempt to an- 
alyze the energyscape whose knowledge is indispensable 
for understanding the dynamics and switching mecha- 
nisms. For this reason, inter alia, calculations of the re- 
versal time of the magnetization of fine single-domain fer- 
romagnetic particles, for instance, initiated by Neel [3,[1], 
and set firmly in the context of the theory of stochastic 
processes by Brown [Oj, ilOj, [2^] , have invariably proceeded 
by ignoring all kind of interactions. Thus, the only terms 
which are taken into account in these calculations are 
the internal magneto-crystalline anisotropy of the parti- 
cle, the random field due to thermal fluctuations, and 
the Zeeman term [ll|. However, before attacking the 
general problem of the effect of interactions and that of 
the surface on the dynamical properties, one has to gain 
a sufficient understanding of the static properties. 

Yet before tackling the effect of surface anisotropy, one 
has to clearly define, understand, and when possible eval- 
uate in a separate way, finite-size effects (artificially) in- 
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duced in theoretical calculations by considering a sys- 
tem of finite size, and "true" boundary effects (naturally) 
induced by crystal field symmetry breaking on the sur- 
face. Finite-size effects in magnetic systems have been 
of much interest for decades now, a.nd have been exten- 
sively studied by many authors ^ IH, IH, [13, S Hi] , 
just to cite a few. One of the difficulties, inherent to sys- 
tems of "round" (spherical or ellipsoidal) geometries, as 
is the case with nanoparticlcs. is the separation of bound- 
ary effects due to defects in coordination on the bound- 
ary and the unavoidable finite-size effects. In hypercubic 
systems, this problem is handled by using periodic and 
free boundary conditions, but this is not possible in other 
topologies, and thus boundary and finite-size effects are 
tangled together. More difhculties come into play when 
the spins on the surface arc attributed some anisotropy. 
While for thin films and layered systems the direction of 
the anisotropy on the surface is clearly determined and its 
intensity estimated ^30.,.3I., i3^. especially in 3d elements, 
for various thicknesses, no definite answer is yet available 
for nanoparticlcs. Nevertheless, there are some estima- 
tions of the surface anisotropy constant [H, [13, [H, [s^ 
obtained in an indirect way by fitting magnetic mea- 
surements and assuming an effective anisotropy constant. 
Unfortunately, today the experimental techniques are not 
so advanced as to allow for probing the crystallographic 
structure of a nanoparticle surface, in the same man- 
ner as is done for 2d magnetism, and hence no infor- 
mation is yet available as to the precise nature of the 
anisotropy on the surface. As such, theorists have no 
definite answer as to what is the most adequate model 
for describing surface anisotropy in nanoparticles. To 
the best of our knowledge there are mainly two models 
with positive or negative constants: The transverse sur- 
face anisotropy (TSA) [2^ model where a surface spin 
is attributed a single-site anisotropy axis that is normal 
to the surface, and Noel's model (NSA) [13, [H, [H, . 
One may also find in the literature a model with random 
anisotropy. Noel's model is more physically plausible be- 
cause the anisotropy at a site occurs only if the latter 
presents defects in its environment, e.g., if it lacks some 
of its neighbors. 

When a nanoparticle is embedded in an assembly dis- 
persed in a non-magnetic matrix (of polymer, silica, etc.), 
still more interesting challenging physical phenomena are 
observed, such as the disappearance of the maximum of 
the temperature at the peak of the zero-field-cooled mag- 
netization of an interacting assembly, the magnetization 
enhancement at low temperature and high field in di- 
lute maghemite particles [l9l . [4l| , and the appearance of 
spin-glass states at low temperature in concentrated as- 
semblies [T9| . In addition to the volume and anisotropy 
easy axes distributions, particle assemblies also introduce 
a new physical parameter, namely the inter-particle in- 
teractions, as discussed earlier, and in particular DDI. 
These interactions have been studied for many decades 
in different areas of physics, and in particular in mag- 
netism [42| . There is an overwhelming literature investi- 



gating the influence of DDI on both static and dynamic 
properties of nanoparticles, and it would be almost im- 
possible to provide a fair account of all of them. How- 
ever, we do not know of theoretical work that deals with 
DDI between particles taking at the same time account 
of their internal structure and thereby surface effects. 
In Ref. [43I . the Landau-Lifshitz thermodynamic pertur- 
bation theory [4^ is used to tackle the case of weakly 
dipolar-interacting monodisperse assemblies of magnetic 
moments, in the OSP approximation, with uniformly or 
randomly distributed anisotropy axes. The authors stud- 
ied the influence of DDI on the susceptibility and spe- 
cific heat of the assembly. On the other hand, Monte 
Carlo technique has been used to study the effect of DDI 
(of arbitrary intensity) on the temperature behavior of 
the coercive field and remanent magnetization in a uni- 
form assembly of magnetic moments (45| . The situation 
involving both surface effects and DDI has never been 
considered so far because of its tremendous complexities. 
Needless to say that, already at the static level, no exact 
analytical treatment of any kind is ever possible even in 
the OSP approximation, i.e., ignoring the internal struc- 
ture of the particles. 

This review is organized as follows: In section II, which 
deals with the intrinsic features of nanoparticles, we em- 
phasize the difference between finite-size, boundary and 
surface effects. Then we introduce our many-spin ap- 
proach and discuss and compare two surface anisotropy 
models (Transverse and Ncel). Next, we present and in- 
terpret some of the main results on the static behavior of 
many-spin particles, and in particular on the role of sur- 
face anisotropy. We end this section with a discussion of 
the crossover between the OSP and MSP regimes. Sec- 
tion II deals with the collective behavior of nanoparticle 
assemblies with special emphasis on DDI: i) we summa- 
rize our work on the effect of anisotropy and weak DDI 
on the magnetization of a polydisperse assembly of mag- 
netic moments, i.e., OSP particles. We provide the cor- 
responding approximate analytical expressions and com- 
pare them with Monte Carlo simulations, ii) We investi- 
gate the effect of DDI on the ZFC magnetization, and in 
particular on the temperature at its maximum as a func- 
tion of the applied magnetic field, iii) We end this section 
with a very short discussion, together with a preliminary 
result, of assemblies of non-interacting MSP particles. 



II. INTRINSIC EFFECTS: SURFACE 
ANISOTROPY 



As discussed in the introduction, when dealing with 
fine magnetic systems, such as a nanomagnct, one should 
distinguish, at least from a theoretical point of view, be- 
tween finite-size, boundary, and surface effects. 
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A. Finite-size versus boundary effects 

For instance, for a simple cubic (sc) lattice (see Fig. [1] 
left) with periodic boundary conditions (pbc), there is 
only one environment (crystal field) with coordination 
number z = 6. In this case, the temperature behav- 
ior of the magnetization is marked by the well-known 
M ~ tail in the critical region, where M is the 

total number of spins in the particle (Fig. [Tl right). In 
the case of more realistic free boundary conditions (fbc) , 
a cube with sc structure shows four different environ- 
ments with z = 3,4,5,6 (see Fig. [5] left). In this case, 
in addition to the previous finite-size effects, one is faced 
with boundary effects. These induce stronger fluctua- 
tions that suppress the magnetization of the system (see 
Fig. [Alright). Considering both cases of pbc and fbc al- 
lows for a separate estimation of the related effects. Now, 
if the boundary of a system with fbc is endowed with 
a surface anisotropy, which is indeed a consequence of 
boundary defects, we may then speak of surface effects, 
in addition to the finite-size and boundary effects. 

For both pbc and fbc, we have shown [i^ that the 
magnetization can be written in a simple form. At low 
temperatures and zero field M deviates from 1, its satu- 
ration value, according to 



(1) 



Therefore, Fig. [3] shows that the coefficient in the 
linear-^ term in Eq. ([T]) is smaller than in the bulk for the 

pbc system and greater for the fbc system. Consequently, 
this means that the boundary effects suppress the mag- 
netization while finite-size effects enhance it. Moreover, 
Fig. [3] also shows that boundary effects render a larger 
contribution than finite-size effects, making the net mag- 
netization below that of the bulk. 



B. Surface effects : MSP approach 

Surface effects are due to the breaking of crystal-field 
symmetry at the boundary of the nanoparticle, and in 
reality they cannot be separated from boundary effects. 
In order to study such effects one has to resort to mi- 
croscopic theories, unlike the macroscopic SW model 
B 3 J capable of distinguishing between different atomic 
environments and taking account of physical parame- 
ters such as single-site surface anisotropy, exchange and 
DDI, in addition of course to the magneto-crystalline 
anisotropy in the core and magnetic field. Unfortunately, 
this leads to difficult many-body problems which can only 
be efficiently dealt with, in general, using numerical ap- 
proaches. 

In this section we will summarize some of the main 
results obtained so far. 



where 



-Ak' 



(2) 



and for a three-dimensional [d = 3) sc lattice Ak = 
(cos kx + cos ky + cos fc^ ) / d. It is important to note that 
Wn in ^ for pbc and fbc differ onl y by the definition of 
the discrete wave vectors, since [46l l47l| 



ka 



l-KriajN-, pbc 
una/N, fbc 



0,1,...,7V-1 (3) 



where a = x,y, z. This subtle difference is responsible for 
much stronger thermal fluctuations in the fbc model due 
to boundary effects. The difference between the sums 
Wn (with finite N, which is the side of the cube) and 
the so-called Watson's integral (for bulk) 



W = 



d^k 



1 



(27r)3 1 - A, 



with W = 1.51639 for d = 3, describes the finite-size 
effects for thepbc case [U, and boundary effects in 
the fbc case [401. Indeed, we have [491 



1. Model and notation 

We consider a ferromagnetic particle of M spins cut 
from a cube of side N (i.e., — 1 atomic spacings) with 
a given lattice structure. Due to the underlying (dis- 
crete) lattice structure, the particle thus obtained is not 
a sphere with smooth boundary because its outer shell 
presents apices, steps, and facets, resulting in many sites 
with different coordination numbers. 

The model Hamiltonian we use is the (classical) 
anisotropic Dirac-Hcisenberg model [H, 



1~L — ^ ^ Jij^i 



(6) 



where is the unit spin vector on site i, H the uniform 
magnetic field, M the total number of spins (core and sur- 
(4) face), and Jij{= J > 0) the nearest-neighbor ferromag- 



netic exchange coupling, 
anisotropy energy 



TLan is the uniaxial single-site 



(7) 



Wn-W 

w 



0.90 

91n(1.17iV) 
2'kWN ' 



pbc 



fbc 



with easy axis and constant Ki > 0. If the spin at site 
(^5^ i is in the core, the anisotropy axis is taken along the 
reference z axis and Ki = Kc- For surface spins, this axis 
is along the radial (i.e., transverse to the cluster surface) 
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FIG. 1: Cubic system with pbc and thermal behavior of the magnetization for two particle sizes. 
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FIG. 2: Cubic system with fbc and thermal behavior of the magnetization for pbc and fbc systems. 



direction and Ki = Kg. In this case, the model in ([7]) for 
surface spins is called the TSA model introduced earlier. 
We use the more general model of transverse direction 
given by the gradient [the vector perpendicular to the 
isotimic surface ■0 = constant defining the shape of the 
particle, e.g. a sphere or an ellipsoid]. In the case of a 
spherical particle the transverse and radial directions co- 
incide, whereas for another geometry such as an ellipsoid 
they do not. 

A more physically appealing microsco pic model of sur- 
face anisotropy was introduced by Neel [37[ with 



n 



Nccl 



(8) 



where Zi is the coordination number of site i and = 
Vij jvij is the unit vector connecting the site i to its near- 
est neighbors j = 1, . . . , z^. This model is more realistic 
since the anisotropy at a given site occurs only when the 
latter loses some of its neighbors, i.e., when it is located 
on the boundary. The model in ([8]) is referred to as the 
Neel surface anisotropy (NS A) model [4§| . 

Qualitatively, the NSA model is not quite different 
from the TSA model. For example, consider a site i sit- 
ting on a, [100] facet, e.g. in the upper most plane normal 
to the z axis. It has 4 neighbors on that facet and one 
below it along the z axis. From ([5]), the corresponding 
energy reads [50| 



n. 



NSA 



2s 



I 



where we have used || s.; ||= 1. This implies that if Kg > the easy direction is along ite^, i.e., normal to the facet. 



6 



3- 










0.2- 


fbc 




9ln(1.17A0/(27rAW) 




0.1 - 










0.0- 


















-0.1 - 


pbc^ 




A}^'^'^' = -0.90W 




-0.2- 










-0.3- 











0.0 0.1 0.2 0.3 j/^ 



FIG. 3; Lattice sums Wn for cubic systems with free and 
periodic boundary conditions. W = 1.51639 is the bulk value 
for the sc lattice. 

and if Ks < the facet becomes an easy plane. There- 
fore, upon dropping the irrelevant constant, we rewrite 
the above energy as 

which is the same as the TSA in Eq. ([7]) for the site con- 
sidered. More generally, averaging the NSA over a sur- 
face perpendicular to the direction n leads to [see Eqs. (6, 
7) in Ref. 

nf^^ = -Ks {WAsl + \ny\sl + \nAsl) , (9) 

thus favoring among {x^y, z) the direction closest to the 
surface normal. This explains the similarity between the 
results obtained with TSA and NSA. Using the compo- 
nents |rict|,Q; — x,y,z, the atomic surface density reads 

m 

/(n) =max{|n^|,|ny|,|n^|}, (10) 
and thereby we have the particular cases 

For comparison, taking account of the atomic surface 
density, the TSA model is described by 

= -/v. (n • s,;) V(n) (11) 

and the whole effect comes from the atomic surface den- 
sity. Quantitatively, in the NSA model the effect is big- 
ger, since for the surface cut perpendicular to the grand 
diagonal of the cubic lattice the anisotropy completely 
disappears, whereas in (jlip the surface anisotropy is only 
reduced by the factor 1 / VS. 



2. Results and discussion 

The static properties of magnetic nanoparticles as 
many-spin systems have been studied in Refs. 
113, I49I . [51I using the classical Heisenberg model in- 
volving the atomic spin, as the elementary building block, 
with continuous degrees of freedom. In these works we 
investigated the field, temperature, and spatial behavior 
of the net as well as the local magnetization. We adapted 
and applied the existing theoretical approaches, such as 
spin-wave theory, spherical model, Monte Carlo simula- 
tion, and numerical techniques to nano-scaled magnetic 
systems. In particular, we studied the magnetic struc- 
ture and found new switching mechanisms of the par- 
ticle endowed with surface anisotropy. It was possible 
to estimate and compare finite-size, boundary, and sur- 
face effects, and to show that the latter two effects are of 
long-range and progressive, and that surface effects in- 
duce non-saturation of the magnetization, as observed in 
experiments [l^[l^, [l^. In the particular (and typical) 
situation where the exchange interaction is much stronger 
than the anisotropy energy we found that the magneti- 
zation can be represented as the spatially homogeneous 
"global" magnetization plus a small inhomogeneous con- 
tribution that we calculated analytically and numerically 
pgl . Isoj . The latter is induced by surface anisotropy and 
it is maximal near the surface but can extend deeply into 
the body of the particle. It describes the adjustment 
of the magnetization to the conditions at the surface by 
minimizing the total energy with fixed direction of the 
global magnetization. As a result we obtained the ef- 
fective particle's energy that depends on the orientation 
of its net magnetization and arises because of surface 
anisotropy. This contribution is of second order in the 
surface anisotropy and it adds to other terms, such as the 
bulk anisotropy and the first-order contribution from the 
surface anisotropy, which disappear in samples of cubic 
or spherical shape. These contributions to the energy of a 
magnetic nanoparticle are crucial to its dynamical behav- 
ior, in particular, in the ferromagnetic resonance (FMR). 
Accurately taking all of them into account should make 
it possible to determine the bulk and surface anisotropics 
from the experimental data. An interesting problem is 
the dynamical aspect of the magnetization adjustment 
mentioned above. If the anisotropics are much smaller 
than exchange interaction, the exchange-driven adjust- 
ment is much faster than the global precession of the 
magnetization induced by the anisotropy. Then, these 
adjustment modes behave adiabatically at low frequen- 
cies and the effective OSP energy is a good approxima- 
tion. For materials with a very strong surface anisotropy 
such separation of dynamical scales is no longer valid, 
and the dynamics of such nanoparticles becomes an es- 
sentially many-body process. FMR experiments on mag- 
netic nanoparticles should allow to estimate the values 
of the surface anisotropy and detect different regimes of 
their dynamical behavior. 

Let us now review with data plots a few typical re- 
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FIG. 4: Hysteresis for ip = 11/4:, M — 360 and different values 
of surface anisotropy constant ks = Ks/Kc, h = H/2Kc- 



suits obtained for an MSP particle with core and sur- 
face anisotropics, together with intra-exchange inter- 
actions and magnetic field. At zero temperature, in 
Refs . [T3I [l^. [51I we studied the magnetic state and switch- 
ing mechanisms of model particles of spherical shape, sc 
lattice, uniaxial anisotropy in the core and both TSA and 
NSA, by varying the size, the exchange couplings, surface 
anisotropy intensity, and the field direction. The main re- 
sult of this work is that varying the surface anisotropy one 
observes a crossover from i) the coherent-reversal regime, 
as obtained by representing the particle as a macrospin 
according to the models of SW (for statics) and NB (for 
dynamics), into ii) the incoherent-reversal regime with 
cluster-wise switching. In the latter regime, the parti- 
cle exhibits (due to strong surface anisotropy) new fea- 
tures that are reminiscent of a many-spin system and 
which cannot be described by a macroscopic approach. 
Fig. [3] shows hysteresis loops for a spherical particle of 
Af = 360 spins with uniaxial anisotropy (with constant 
Kc) in the core and TSA with constant kg = Kg/J. 
The field is applied at an angle i/) = 7r/4 with respect 
to the core easy axis. These results show that as the sur- 
face anisotropy constant increases, one starts to observe 
steps in the hysteresis loop corresponding to the succes- 
sive switching of spin clusters, defined as the groups of 
spins whose anisotropy axes make the same angle with 
the field direction. Such subtle (intrinsic) features could 
be checked by direct measurements of the hysteresis loop 
of an individual nanoparticle, were they to become possi- 
ble. At present, available measurements of nanoparticle 
assemblies do not allow for such a check. 

Fig. O shows the 2D limit-of-metastability curve (or 
SW astroid) for the particle as in Fig. 21 On the left 



panel we see that for small surface anisotropy this curve 
perfectly scales with the SW model (in full line) with 
the scaling factor Nc/M, where Nc is the number of core 
spins, i.e. those spins with full coordination number. On 
the other hand, the right panel shows that as the surface 
anisotropy increases this curve depresses in the longitu- 
dinal direction, thus deviating from the SW curve. It 
is clear that as Ks increases the spin non-coUinearities 
arc enhanced and the SW model is no longer applicable. 
In Ref. m we showed that the spin disorder caused by 
boundary defects is of long range and propagates to the 
center of the particle, which is at variance with the so- 
called core-shell model, as far as the magnetic properties 
arc concerned. We have also studied minor hysteresis 
loops as shown in Fig. [51 Their appearance is an indica- 
tion of the existence of several local minima induced by 
surface anisotropy. Indeed, we found that the bigger is 
the particle, and thereby the smaller is the surface contri- 
bution, and the narrower is the minor loop, which implies 
that the non-coUinearities caused by surface anisotropy 
are indeed responsible for the minor loop. 

At finite temperature, one can study equilibrium ther- 
modynamic properties by Monte Carlo simulations (48j 
or b y th e adapted spherical model [l^ or spin-wave the- 
ory [431. For a box-shaped particle with pbc or fbc and 
without surface anisotropy, we computed the magnetiza- 
tion as a function of field and temperature both analyt- 
ically and numerically. We used respectively the modi- 
fied spin-wave theory and Monte Carlo simulation using 
an augmented Metropolis algorithm that takes accounts 
of the superparamagnetic global rotation of the whole 
bunch of spins in the particle. For instance. Fig. [71 shows 
the magnetization of a many-spin particle as a function 
of magnetic field at different temperatures. These re- 
sults obtained by Monte Carlo simulations simply show 
that because of surface effects the magnetization docs 
not saturate even at relatively high fields, i.e. h = 0.2 
which corresponds to ~ 16 Tesla for cobalt particles 
even at very low temperature, i.e. t = 0.01 (T ^ 1 K). 
In Ref. [13, it was shown that the long-range order of 
the spin disorder at the particle's boundary, mentioned 
above, is accentuated by thermal fluctuations. 

The study of dynamic switching of the magnetization 
triggered by thermal fluctuations is more involved and 
one has to resort to stochastic approaches based on the 
Landau-Lifshitz-Langevin equation, as a dynamic equa- 
tion for the magnetization, or the Fokker-Planck equa- 
tion, as a dynamic equation for the probability distribu- 
tion of the magnetization orientations. It is also possible 
to use the so-called kinetic Monte Carlo technique. On 
the other hand, the approach based on the solution of the 
Fokker-Planck equation by the matrix-continued fraction 
method [s^l is limited to a small number of degrees of 
freedom, or equivalently, a small number of spins inside 
the particle. The main reason is that with a many-spin 
particle one is faced with complex many-body aspects 
with the inherent difficulties related with analysing the 
energyscape (location of the minima, maxima, and sad- 
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FIG. 5: Astroid of a spherical particle of A/" = 360 spins. Left: for ks — Ks/Ka ~ 1 and various values of the surface-to- volume 
ratio Nst = Ns/M. Right: various values of surface anisotropy constant ks = Ks/Kc- The full dark line is the SW astroid 
scaled with Nc/Af, where Nc is the number of spins in the core; the dotted line is only a guide to the eye. 
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FIG. 6: Minor hysteresis loops for the particle diameters D = 
9, 10, 11, 13, which corresponds to the sizes indicated in the 
legend. Kc/J = 0.01, Ks/ J = 0.3. The return field is h ^ 
-0.4 {h = H/Hk, Hk = 2K,). 



FIG. 7: Magnetization as a function of field with variable 
temperature {t = knT/J) for a particle of 515 spins with core 
uniaxial anisotropy (K^/J ~ 0.0024) and transverse surface 
anisotropy {Ks/J — 0.4). 



die points of the energy potential). This analysis is, how- 
ever, unavoidable since it is a crucial step in the calcu- 
lation of the relaxation time and thereby in the study of 
the magnetization stability against thermally-activated 
reversal. On the other hand, the Langevin approach is 
quite versatile as it can be applied to a general system, 
but then one has to average over many stochastic trajec- 
tories generated with a random Gaussian field. This ren- 
ders this approach rather time consuming. For an OSP, 
the limit-of-mctastability curve has been obtained at dif- 
ferent temperatures in Ref. [H. Such thermal effects had 
also been studied in Ref. [13 by constructing a free energy 
that interpolates between the low-temperature Landau- 
Lifshitz energy and the high-temperature Landau energy. 
It was found that there is a crossover between the low- 
temperature regime where the magnetization switching 
follows the SW model, i.e., through coherent rotation, 
and the high-tcmpcrature regime where the magnetiza- 



tion switches by varying its magnitude. In both works, it 
is shown that the thermal fluctuations induces a shrink- 
ing of the limit-of-metastability curve, which was ob- 
served by ^-SQUID measurements 0. In the case of 
a many-spin particle Landau-Lifshitz-Langevin calcula- 
tions may be performed in order to investigate the ef- 
fects of both temperature and surface anisotropy on the 
hysteresis loop and limit-of-metastability curve. 

In Ref. we investigated, in a somewhat phenomcno- 
logical way, the effect of varying the ratio of the exchange 
coupling at the surface to that in the core in maghemite 
particles. The idea was to model the influence of the ma- 
trix in which the particles are embedded. The exchange 
coupling at the surface is meant to be an effective cou- 
pling resulting from the various chemical and physical 
interactions between the matrix and the particles. Us- 
ing the Monte Carlo technique we computed the thermal 
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variation of the surface magnetisation as a function of 
the (reduced) temperature for a ratio of the exchange 
couphng at the surface to that in the core Jsc = 0.5, 1, 2 
and for a surface contribution of N^t = 40%. We found 
that the surface critical region is shifted to higher tem- 
peratures upon increasing J^, and only when Jg = 2Jc 
do both the core and surface magnetic phase transitions 
occur in the same temperature range. The results also 
showed that the weaker is the exchange interactions on 
the surface the lower is the magnetisation of the latter. 
This result remains the same upon decreasing the outer 
shell thickness. In this case the number of spins having 
smaller coordination numbers than in the core, and hence 
weaker effective exchange energy, i.e., those spins on the 
outer shell of the particle, is very small as compared to 
the rest of spins in the particle. Moreover, wc may think 
that the spins on the outer shell follow the (strong) effec- 
tive field created by the other (inner) spins constituting 
a relatively ferrimagnetically ordered core. 

C. MSP versus effective OSP 

With the desire to avoid the difficulties mentioned 
above, one may address the question as to whether there 
exist some cases in which the full-fledged theory that 
has been developed for the OSP approach [see [54I and 
references therein] can still be used to describe the dy- 
namics of an MSP. However, avoiding somehow the spin 
non-coUincaritics induced by surface/interface anisotropy 
means that some price has to be paid. In Rcf. [1^ it was 
shown that when the surface anisotropy is much smaller 
than the exchange interaction, and in the absence of core 
anisotropy, the surface anisotropy contribution to the 
particle's energy is of 4"^-order in the net magnetization 
components and 2'''^-order in the surface anisotropy con- 
stant. This means that the energy of an MSP, cut from 
a lattice of cubic symmetry, with relatively weak surface 
anisotropy can be modeled by that of an OSP whose effec- 
tive energy contains an additional cubic-anisotropy po- 
tential. An analytical expression was given for the effec- 
tive constant Kcs of the surface-induced cubic-anisotropy 
term, when the core anisotropy is absent, that is 

= (12) 

where JV, Ks,z,J are respectively the number of atoms, 
the surface anisotropy constant (transverse or Neel), the 
coordination number, and the exchange coupling of the 
many-spin particle, k is a surface integral that depends 
on the underlying lattice, the shape, and the size of the 
particle and also on the surface-anisotropy model. For 
a spherical particle (of ~ 1500 spins) cut from a simple 
cubic lattice with Neel's surface anisotropy, k ~ 0.53465. 
In the presence of core anisotropy the effective energy 
contains terms with coefficients that are products of the 
core and surface anisotropy constants. However, for small 
surface anisotropy the effective energy again provides a 



good approximation of the many-spin particle, as has 
been shown in Ref. [lOl. In this case the energy of a 
many-spin particle with uniaxial anisotropy in the core 
and TSA or NSA on the surface is modeled by that of a 
one-spin particle with the net magnetization m in a po- 
tential containing a uniaxial and cubic anisotropy terms, 
i.e., up to a constant, wc have 

f off = - A'uni ml + Kcuh{ml +m^ + ml). (13) 

The results are shown in Fig. [8l 

Before wc discuss the results let us briefly explain the 
method we used to obtain them. Because we are dealing 
with an MSP, the energyscape cannot be represented in 
terms of the coordinates of all spins. Instead, we may 
represent it in terms of the coordinates of the particle's 
net magnetization. For this purpose, we fix the global or 
net magnetization, m, of the particle in a desired direc- 
tion mo (I^Tiol = 1) by using the energy function with a 
Lagrange multiplier A [i^ : 

V s 

J=-^H-AfX-im~mo), m = (14) 

To minimize J^, we solve the evolution equations 

s, = - [si X [si X F^]] , Fi = -dJ^/ds, 

X = dT/d\^-JV{m-mo), (15) 

starting from Si ~ mo = m for all i = 1, . . . ,J\f and 
A = 0, until a stationary state is reached. In this state 
m = mo and [s^ x F^] = 0, i.e., the torque due to the 
term JVX- (m — mo) in compensates for the torque 
acting to rotate the global magnetization towards the 
minimum-energy directions [see discussion in Ref. [49j]. 
The orientation of the net magnetization is then given 
either in Cartesian coordinates {uix, niy, rriz) or in spher- 
ical coordinates {On, ^n)- 

The first panel of Fig. M {ks = K^/J = 0.1) shows 
that for very small kg the energyscape of an MSP is 
well recovered by the effective energy in Eq. p^ . As 
kg increases [see middle panel, kg = 0.3], some devia- 
tions start to be seen, and for relatively large values of 
kg a fit with Eq. p3)) is no longer possible. In fact, in 
this regime strong deviations from coUincarity develop, 
especially near maxima and saddle points. In fact, in 
this case the Lagrange-parameter method of Eqs. (HH 
[TS)) fails because the magnetic state of an MSP can no 
longer be represented by a net magnetization. These re- 
sults imply that the effect of spin non-coUinearities on the 
energy is to split the minimum at 0„ = 0, defined by the 
core uniaxial anisotropy, into four minima at 0„ ^ 28° 
and (pn = 0, ±7r/2, tt, reminiscent of cubic anisotropy [see 
Fig. [9] . These minima arc connected by saddle points at 
(fin = ±7r/4 and ±37r/4 and the point at 6'„ = be- 
comes a small local maximum. The four minima exist 
over a finite range of the applied field, although their 
positions change continuously as a function of the field 
[5(| . Fig.[Tn]is a plot of the 2D energyscape for a spher- 
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FIG. 8: (Color online) 2D energyscape for the spherical particle of TSA with varying constant kg — Kg/ J. The symbols are 
results of the numerical calculations for the MSP, and the full lines are fits using Eq. (|13p . (fn = and h = 0. 




FIG. 9: Plot of cos (pn, ^ sin 1^3,1, E) for the same param- 
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FIG. 10: (Color online) Energy as a function of the polar 
angle On for a spherical particle with uniaxial anisotropy in 
the core and Neel anisotropy on the surface, ipn = 0, 7r/4 and 
ft = 0. 



ical particle with uniaxial anisotropy in the core, as be- 
fore, but now with NSA on the surface. It is clear that 
the cubic-anisotropy features are also seen in the case of 
NSA, namely that i) the energy minima are not along 
the directions (0„ = 0,7r) of the core easy axis; these 
directions having become local maxima, as discussed in 
the case of TSA, and ii) there is a clear dependence on 
the azimuthal angle ipn- In addition, we would like to 



mention that more extensive calculations [54| have shown 
that similar features are also observed for other crystal 
structures (fee, bcc, etc.). 

These results agree and complement those of Ref. [Tsl 
and [51I. where it was shown that for both the TSA and 
NSA models there exists a (different) critical value of the 
surface anisotropy constant that separates i) the OSP 
SW regime of coherent switching and ii) the MSP regime 
where the strong spin non-coUinearities invalidate the 
coherent mechanism, and the particle can no longer be 
modeled by an effective OSP. Obviously, for very small 
surface anisotropy the cubic contribution becomes neg- 
ligible [see the first panel in Fig. [H fc^ = 0.1] and thus 
the SW OSP model provides a good approximation to 
the many-spin particle. Accordingly, some experimental 
macroscopic estimations of the surface anisotropy con- 
stant yield, e.g^ for cobalt Kg/ J — 0.1 [IHl, for iron 
Ks/J ~ 0.06 [5y|, and for maghemite particles Kg/ J ~ 
0.04 [s^]. However, one should not forget that this ef- 
fective constant depends on the particle's size, among 
other parameters such as the material composition, and 
for, e.g. a diameter of 2 nm we may expect higher 
anisotropies. 

Using the effective energy we have also studied [111 
the ferromagnetic resonance in systems with competing 
uniaxial and cubic anisotropies and built a model that ap- 
plies to i) magnetic materials with both uniaxial and cu- 
bic anisotropies, and ii) magnetic nanoparticlcs with core 
and effective surface anisotropies. We numerically com- 
puted the resonance frequency as a function of the field 
and the resonance field as a function of the direction of 
the applied field for an arbitrary ratio of cubic-to- uniaxial 
anisotropy. We also provided some analytical approxi- 
mate expressions in the case of weak cubic anisotropy 
and proposed a method that uses these expressions for 
estimating the uniaxial and cubic anisotropy constants, 
and for determining the relative orientation of the cubic 
anisotropy axes with respect to the crystal principle axes. 
This method is applicable to the analysis of experimental 
data of resonance type measurements. 

It would very interesting to investigate the crossover 
between the OSP and MSP regimes, in terms of the 
particle's size and environment and material's proper- 
ties. In addition, by comparing the experimental results 
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with various models (e.g., TSA and NSA) for surface 
anisotropy, one should be able to determine the most 
appropriate model. 



III. COLLECTIVE EFFECTS 

While measurements of isolated single particles by the 
technique of /i-SQUID @| have been so successful and 
have provided us with invaluable information on their 
intrinsic static and dynamic properties, they still suffer 
from some limitations. Indeed, apart from the limited 
range of temperature and applied field, the magnetiza- 
tion value is still inaccessible and only the field at which 
it changes sign is detectable. This prevents us from ob- 
taining the magnetization as a function of the field, tem- 
perature and time, and in particular no hysteresis loop is 
available. On the contrary, measurements of particle as- 
semblies are possible in larger ranges of applied field and 
temperature. However, the analysis of the correspond- 
ing results is made rather difficult and subtle because of 
the volume and easy-axis orientation distributions, and 
more importantly, the DDI, whose intensity depends on 
the sample concentration. At present, most of the stud- 
ies on magnetic nanoparticlcs are conducted on assem- 
blies. This is mainly due to the above mentioned rea- 
sons and also because today one may understand better 
how assemblies of nanoparticles can be used in techno- 
logical applications than one does for an isolated par- 
ticle. These applications require ever denser assemblies 
and thus smaller particles. However, this brings a new 
dilemma because small particles become superparamag- 
netic at even low temperature. Moreover, high density 
entails strong DDI among the particles, and in techno- 
logical applications such as magnetic recording, this is 
an issue of special importance because DDI have been 
widely recognized as being responsible for the deterio- 
ration of the signal-to-noise ratio [see e.g.. Refs. [s^, [g^ 
and references therein]. As such, an optimum material 
[with appropriate anisotropy and other physical parame- 
ters] has still to be devised. On the other hand, the study 
of nanoparticle assemblies brings new headaches to the- 
orists, at least, since they are faced with tremendous dif- 
ficulties related with DDI between particles and their in- 
terplay with surface effects intrinsic to the particles. Nev- 
ertheless, the problem of DDI in nanoparticle assemblies 
has triggered much interest due to many new phenomena 
that emerge from the collective behavior of the particles, 
notably the so-called spin-glass state at low temperature 
in concentrated assemblies [see Refs. [H, [6l|, |62|, and 
many references therein] , and also because these interac- 
tions have always constituted a challenging issue in many 
areas of physics. 

Today, there arises the important issue about assem- 
blies of nanoparticles that concerns the understanding 
of the interplay between the intrinsic properties, such as 
those pertaining to surface effects, and extrinsic or collec- 
tive effects stemming from the long-range DDI. Many re- 



search groups have experimentally studied this interplay 
in cobalt and maghemite particle assemblies. Measure- 
ments of the magnetization at high fields per formed on 
the 7-Fe203 nanoparticles [ill, [6^, [see also |l6| for cobalt 
particles] have shown that the magnetization is strongly 
influenced by the particle's size. For instance. Fig. 1 of 
Ref. [i^ shows that i) there is a sudden increase of the 
magnetization as a function of the applied field when the 
temperature reaches 70 K, and the magnetization does 
not saturate at the highest available field, i.e., 5.5 T. ii) 
there is an important increase of the magnetization at 
low temperature, iii) the thermal behavior of the mag- 
netization at 5.5 T is such that the smaller is the mean 
diameter of the particle the faster is the increase of the 
magnetization at very low temperature. 

On the theory side, the situation involving both sur- 
face effects and DDI has never been considered so far 
mainly because of its tremendous complexities, and also 
because one has first to understand these two effects sep- 
arately. Needless to say that, already at the static level, 
no exact analytical treatment of any kind is ever pos- 
sible even in the OSP approximation, i.e., ignoring the 
internal structure of the particles and thereby surface 
effects. Only numerical approaches such as the Monte 
Carlo technique can be of some rescue here. This tech- 
nique has been used in Ref. [i^ to study hysteretic prop- 
erties of monodisperse assemblies of nanoparticles with 
the more realistic Heisenberg spin model, in the OSP ap- 
proximation. In the case of weak interactions, as is ap- 
plicable in dilute assemblies, some approximate analyti- 
cal expressions can be obtained. Indeed, in Ref. 0, the 
Landau-Lifshitz thermodynamic perturbation theory (4^ ] 
is used to tackle the case of weakly dipolar-interacting 
monodisperse assemblies of magnetic moments, i.e., in 
the OSP approximation, with uniformly or randomly dis- 
tributed anisotropy axes. The authors studied the influ- 
ence of DDI on the susceptibility and specific heat of 
the assembly. In Ref. [g^, the same approach was used 
with the objective to study the effect of anisotropy and 
(weak) dipolar interactions on the field and tempera- 
ture behavior of the magnetization of a monodisperse 
and polydispersc assemblies of magnetic moments. We 
studied an assembly of magnetic moments whose magni- 
tudes are distributed according to a lognormal function. 
The anisotropy is taken as uniaxial and either textured 
along some reference axis or randomly distributed. The 
low-field regime magnetization obtained in Ref. [65l . is 
generalized so as to take account of polydispersity and 
DDI. In high fields, the magnetization as a function of 
temperature and field was computed using the steepest- 
descent approximation. In the general range of temper- 
ature, field, and anisotropy, the magnetization of a non- 
interacting assembly was computed exactly by numerical 
integration of the single-moment (free) partition func- 
tion. For interacting assemblies, we used the Monte Carlo 
technique. It was possible to obtain practical (semi) ana- 
lytical formulae for the field and temperature dependence 
of the assembly magnetization that take into account mo- 
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mcnt and easy axes distributions, and weak DDI. It was 
then possible to investigate the effect of anisotropy and 
DDI and to discuss the vahdity of the Langevin law, for 
both textured and random-anisotropy, which is very of- 
ten used in the literature to interpret the magnetization 
measurements on nanoparticle assemblies. 

In the following sections, we review the main results of 
this work and also discuss its extention to include intrin- 
sic effects. 



with I Si I = 1. The magnitude of the magnetic moment 
m.i is defined in terms of the Bohr magneton ^b, i-C, 
rrii ~ riiiiB, and the numbers are either all equal 
for monodisperse assemblies or lognormal distributed for 
polydisperse assemblies. Each magnetic moment is as- 
signed a uniaxial easy axis Bi, and for an assembly these 
axes are either all directed along some reference axis lead- 
ing to a textured assembly, or randomly distributed. 



A. Assembly of OSP particles: effects of anisotropy 
and DDI 

1. Notation and basic equations 

We consider an assembly of magnetic moments — 
rriiSi, i = 1,...,A/' of magnitude mi and direction Si, 



The energy of a magnetic moment ra.i with uniaxial 
anisotropy axis e^, interacting with all the other moments 
via DDI, in the magnetic field H = Heh, reads [after 
multiplying by = -l/fcsT], 



(16) 



j<i 



0<l 



where 



(17) 



is the vector joining the sites i and j and whose magni- 
tude is measured in units of a, a characteristic length on 
the lattice to be specified later, sf'^ is the free particle 
energy, and 



^J'BI< 



(18) 



are the dimensionless energy parameters. It is also con- 
venient to introduce the parameters Xi = xrii , ai = mii . 
Note that ai = KVi/{kBT) is the commonly used no- 
tation for the reduced anisotropy-barrier height of the 
particle i. 

The magnetization per particle of the assembly, or its 
component in the field direction, which is taken here 
along the z axis, is given by 



1 f tP ^ 



In Eq. (dni), (™f) is the statistical thermodynamic aver- 
age over the particle's moment direction s^. 



where Vn = J], d^t = U, d'^s^/2TT, and £ = J2t 



2. Effect of anisotropy and DDI on M{T, H) 



(21) 



{ml) (TOi,e,,cr,x,Cd), 
(19) 



Applying the thermodynamic perturbation theory [44 ] 
and noting that i) since the field is applied along the 
z axis the average of the x, y components vanishes, ii) 
DDI only involve pairs of distinct indices, we obtain the 
following expression for the magnetization of a weakly 
interacting assembly (to first order in ^d) H, 

{S!) ^ {S!),+£,dY.{Sl),A j^^^['^\ (22) 



k=l 



where S,- — riiSi and 



where 'w{ni) is the lognormal distribution of the numbers 
rii with parameters /i, S, 



[3(e,, • ekif - 1] 



dx 



eji ■ Vki ■ eh, 



w{n) 



exp 



1 /Inn — fi 
2 



' kl 



(20) 



— ^3 (^GySy 1) 



(23) 



are the DDI tensors. As was discussed in Ref. and 
confirmed in Ref. [g^, for non spherical systems, the cor- 
rections to the magnetization are largely dominated by 
the first order contribution to the DDI. 

The simple equation ((22|) tells us that the magnetiza- 
tion of a (weakly) interacting particle is given in terms 
of the magnetization (5*^)0 of the free particle and its 
derivatives with respect to the field. {S^)q is the statis- 
tical thermodynamic average in the absence of DDI and 
is thus given by Eq. (PTjl with £i replaced by [see 
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Eq. IHl)]. 



Now, the free particle magnetization can be computed 
either exactly by numerical integration in (|2ip or analyt- 
ically in some limitting cases, such as low and high field. 
Approximate analytical expressions can be obtained in 
the low-field regime by perturbation theory [g^, [sl, [6^ 
and in the high-field regime using the steepest-descent 
approximation [6^ . More precisely, in low fields we have 



Xi — textured 

(24) 



■ 2 

Xi 1 + 2S^2 3 



3 45 



random, 



where [6£ 



{21- 



1 - 



^(f)'/2 + ..., a.«l, 



1(1 + 1) 

4(Ti 



The high-field magnetization reads [6; 

1 , 

I ^ ~ 

\hf 



CT, > 1. 



textured 



1 af 1 
1 E"~2' random. 

Xi ±0 X- 



(25) 



(26) 



The validity of the asymptotic low field and high field 
expressions of the magnetization given in and (pS)) . 
respectively, is checked by comparing them with the ex- 
act numerical calculation of the partition function in 



(PT|) . We find that these asymptotic expressions are good 
enough even at a relatively high temperature (here small 
a), for which the steepest-descent approximation is ex- 
pected to work worst because x becomes small. More 
general expressions were obtained in Ref. [g^ by a differ- 
ent approach. 



Next, we generalized the expressions analogous to 
Eqs. and to a weakly interacting polydispcrse 
assembly in the case of random anisotropy. Accordingly, 
inserting the low and high field expansions ([M]) . (|26p in 
Eq. (|22[1 leads to analytical expressions for the magneti- 
zation of a weakly interacting assembly as a function of 
field, temperature, anisotropy, and the DDI parameter 
^d- Therefore, the assembly (reduced) magnetization per 
particle defined as (s^)^^^ = ^ /J^ YltLi z) i, reads (in 
the case of randomly distributed easy axes) 



3 



3 



^3 + g^dAs 



45 



1 - 



15 



2(a) 
15 



low field 



high field. 



(27) 



where we have defined [see Eq. (fT^l) for notation] l/A^Si;=i ^^'^ introduced the (scaled) asscmbly-on- 
|d = id{nf ,{x) = {n)x,{a) = (n) cr, with (n) = lattice constants 



14 



1 



E 



\a+b 



Z — 1 



1 



E 



(28) 



a.; = 1 + 25j^2 '^it^ '5'i2 being defined in ([25 



In Eq. (j27p . we observe that there are three types of 
contributions: There are pure anisotropy terms, pure 
DDI terms, and mixed terms. It is readily seen that in the 
absence of anisotropy and DDI, i.e., for ai = 0,£,d ^ 0, 
and from Eq. (|25p = 1 , expressions ([??)) simplify back 
to Eqs. fM]) . The presence of mixed terms sim- 

ply reminds us that DDI induce an additional anisotropy 
in the magnetic system and thus re-normalizc the pure 
magneto-crystalline anisotropy. 

In the continuum limit the lattice sum C^^''^^ = 
^ EZ=iM, becomes C(o.o) = 47r(l/3 - A,), for 
a simple cubic lattice, with being the demagnetizing 
factor along the z axis. In the monodisperse case we have 
Q{a,b) ^ ^(0,0)^ rpj_^^ constants C'^'''''\ and thereby the 
corresponding DDI terms in Eq. (j27p . arc shape depen- 
dent, which is no surprise knowing that the long range 
DDI lead to shape dependence of the physical quanti- 
ties, and in particular the magnetization. On the other 
hand, we found that these constants are negative for the 
oblate system, and positive for the prolate, which im- 
plies that the DDI suppress the assembly magnetization 
in the former case and enhance it in the latter. Moreover, 
for cubic systems all the constants C^"'''-' vanish, which 
means that the DDI do not contribute to the magneti- 
zation in this case (see Fig. 4 of Ref. [63 ). and thus the 
deviations from the Langevin behavior are caused only 
by anisotropy. Note that all the scaled constants C^"'''-' 
are almost independent of the assembly mean diameter. 
In particular, this is trivial for C^^'^\ The constants 
^3,^5 contain the parameter and thereby depend 
on anisotropy. From Eq. (|25p we infer that in the ab- 
sence of anisotropy, i.e., Ci = 0, a,; = 1, while for strong 
anisotropy we may approximate Si2 to 1, and hence to 
3, so that ^3 oc (1/7V) J2i "f/ (^>^ and A5 cx C^^.s) in the 
two limits of anisotropy. For Dm = 3,7 nm, ^3 ~ 2,6. 
In the continuum limit A3 tends to exp(35^), where 6 is 
the standard deviation of the distribution A5 ex- 

hibits the same behavior as C^^'^' but with bigger change 
with Dm.. It is well known from other areas of physics 
that the calculation of such high-order moments (or "cu- 
mulants") requires more precision because they present 
more statistical fluctuations with the lattice size. 

In Fig. [11] we plot the Langevin function (full line) and 
the Monte Carlo results (symbols) for the magnetization 
of an interacting assembly of {Af = lOx lOx 5) lognormal- 
distributed moments, with random anisotropy, and for 
three values of the inter-particle distance. The intensity 
of DDI, or cquivalently the value of is varied by vary- 
ing the lattice parameter a entering S,d [see Eq. ([T| 



More precisely, the parameter a is taken as a real num- 
ber times the mean diameter Dm of the assembly, i.e., 
a = k X Dm- Thus, large values of k correspond to an 
isotropically inflated lattice with large distances between 
the magnetic moments, and thereby weaker DDI. These 
results, obtained for an assembly on a simple cubic lat- 
tice, do confirm that DDI suppress the magnetization. 
Indeed, we recall that it was shown by Luttinger and 
Tisza |67[ [see also the more recent work |45| using the 
Monte Carlo technique] that the ground state of a simple 
cubic lattice of dipolcs is antiferromagnetic, while that of 
a face-centered cubic lattice is ferromagnetic. It is also 
seen that these curves deviate from the Langevin law. 
However, we emphasize that the deviations induced by 
DDI are much smaller than those induce d by anisotropy, 
as already discussed earlier [see also Ref. |4j for a related 
discussion of the effect of the system shape on the mag- 
netic susceptibility of a monodisperse assembly]. In the 
inset of Fig. [n] the same results are magnified by plotting 
them in function of C = {{x) /^(i)xlO~^ oc fisH/^fi'^g/a^), 
i.e., the ratio of Zeeman energy to the DDI energy, which 
also makes it possible to distinctively plot the analytical 
expressions ^7} for low fields. In the case of high fields 
only one curve (fc = 2, i.e., relatively strong DDI) is 
presented since for the other values of DDI parameter fc, 
the steepest-descent approximation is only valid for much 
higher values of (. Note also that in function of the pa- 
rameter C the tendency with increasing DDI strength is 
reversed. 



3. Effect of DDI on the ZFC magnetization of an OSP 
assembly 

To analyze and understand the dynamics experiments 
on interacting nanoparticle assemblies, and in particular 
to understand the system's dynamic response such as the 
ac susceptibility and zero-field-cooled (ZFC) and field- 
cooled (FC) magnetization, we need to know how the 
DDI affect the relaxation time of the assembly. While a 
fair understanding of the mechanisms behind the ZFC- 
FC magnetization and ac susceptibility has been achieved 
in the case of non-interacting assemblies, many exper- 
imental results on interacting assemblies remain unex- 
plained, notably the observation, mentioned in the intro- 
duction to section Um of a glass-like state at low tempera- 
ture in strongly interacting assemblies. Obviously, this is 
mainly due to the complexities of the DDI and also of the 
very calculation of the relaxation time itself. Recently, 
Jonsson and Garcia-Palacios (see also [6^) have es- 
tablished an approximate expression for the relaxation 
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Reduced magnetization 




FIG. 11: Reduced magnetization (per particle) of an interacting assembly of A/" = 10 x 10 x 5 lognormal-distributed magnetic 
moments with mean diameter Dm = 7nm and random anisotropy. Monte-Carlo in symbols and analytical expressions (|27p in 
lines. The parameters C, and k are defined in the text. 



rate of a weakly interacting monodispcrse assembly with 
textured or randomly-distributed anisotropy. Then us- 
ing the simple Debey relaxation model they investigated 
the effect of weak DDI on the ac susceptibility and in 
particular the displacement of the maximum of its real 
and imaginary parts. In their explanation of the results 
they emphasized the important role played by damping 
in the relaxation processes in the presence of a transverse 
field, in addition to the effect of the change in the energy 
barriers, that was commonly believed to play the major 
role. The role of a transverse field is played here by the 
transverse component of the dipolar field. 

Experimental results obtained for ferrofluids [tO] and 
later for 7-Fe203 nanoparticles [7l| indicated that for di- 
lute samples (weak DDI), the temperature Tmax at the 
maximum of ZFC magnetization first increases with in- 
creasing field, attains a maximum and then decreases. 
More experiments performed on the 7-FC2O3 particles 
dispersed in polymer [4l| confirmed the previous results 
for dilute samples and showed that, on the contrary, for 
concentrated samples (strong DDI) Tmax was a mono- 
tonic decreasing function of the magnetic field. In Ref.pT^ 
we showed that the bell-like shape of Tmax(^f), in the di- 
lute case, is insensitive to the intrinsic properties of the 
particles, of course in the OSP approximation. Exact nu- 
merical calculations (TSLpT^lTSj of the smallest eigenvalue 
of the Fokker-Planck matrix invariably led to a mono- 
tonic decrease in the blocking temperature, and thereby 
in the temperature Tmax, as a function of the magnetic 
field. Indeed, in Ref. [t^ we showed that the expression 
of the single-particle relaxation time does not seem to 
play a crucial role and that even the (relatively) sim- 
ple Neel-Brown expression for the relaxation time leads 



to a maximum in Ti„ax(-ff) if one considers an assem- 
bly of particles whose magnetization, formulated through 
the Gittleman-Abeles-Bozowski model [t^ , has a super- 
paramagnetic contribution that is a non-linear function 
(such as Langevin's) of the magnetic field. The magneto- 
crystalline anisotropy and the volume-distribution width 
too have strong influence. The issue of the effect of DDI 
on Tmax(-ff)j namely the disappearance of the maximum 
when the intensity of interactions increases, was left open 
in Ref. [t^- In a recent work [tJ we revisited this is- 
sue after generalizing the work of Ref. [6^ to include the 
magnetic field and magnetic-moment distribution (poly- 
disperse assembly) of nanoparticles in the OSP approx- 
imation. We then investigated the effect of (weak) DDI 
on the ZFC magnetization, in particular on Tmax(-f^)- 
Below, we briefly outline the main results [77j . 

The idea is to introduce the effective field ^ composed 
of the external magnetic field and the dipolar field (re- 
member that Vii = and see notation in section [III A 1[) . 

Ci = XjG/i + ^f'P = Xieji + ^ I?ijSj. (29) 

i 

Then the relaxation rate of the assembly is computed 
as the average of the relaxation rates of each magnetic 
moment that experiences this effective field using per- 
turbation theory assuming that | |^ 1- In Ref. [tI 
an estimation of DDI is given for two samples of Cobalt 
nanoparticles which indicates that the DDI field is of the 
order of 300 Oe, which in reduced units, obtained after di- 
viding by the corresponding anisotropy field of the order 
of 0.3 T, is ~ 4 X 10-3 - 10-2. This is of course very 
small. On the other hand, the field at which T,„ax(-ff) 
has a maximum is circa 100 Oe [see Ref. [t^ Fig. 1, for 
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maghcmitc particles] , corresponding to a reduced field of 
3 X 10^^. This indeed suggests that in typical samples 
the above condition on the effective field C is satisfied. 

The relaxation rate of a weakly DDI-interacting nano- 
magnet obtained in Ref. [6^ is written as 



.(0) 



(30) 



where r^'''' is the relaxation rate of the nanomagnet at 
site i in the absence of the DDI field and is given by 



where fiV) is the (normalized) lognormal volume distri- 
bution with parameters Vq and d 



VV = f(y)dV 



1 



exp 



2^2 



dV 
'V' 



Xv is calculated by assuming a step function for the mag- 
netic field, i.e. H = for t < and H = Hq = const for 
t > 0. Then, the contribution to the magnetization from 
particles of volume V is given by 



.(0) 



1/2 - 

=cr/ e 



(31) 



Mvit)^VHo\xo-{xo-Xi)e-'/' 



(35) 



with Ts = {-fHa) ^ which evaluates to 7 ~ 1.76 x 10" (T. 
s)~^ for cobalt particles with the anisotropy field Ha ~ 
0.3 T; 7 ~ 1.76x 10" (T. s)-^ In the high-energy barrier 
approximation, a ^ 1, h = x/2a <^ 1, which we will be 
using later on, the function F reads [73] 



F{a^) ~ 1 



5 1 



(32) 



where A is the Landau-Lifshitz damping parameter. The 
field ^ in Eq. ([SO]) was then replaced by the average over 
the assembly index and spin orientations of the dipolar 
field. 

In Ref. [73, we replace the field by the effective field 
Ci in Eq. ([29]) and then compute the statistical average 
over the spin orientation and the (random) anisotropy 
easy-axis distribution and obtain the double averages 



1 2 



3 ' 



fnf (2i?,,) = 2(c2, 



where 



is a lattice sum. 

Therefore, the relaxation rate of a weakly interacting 
particle containing n.i Bohr magnetons, embedded in a 
polydisperse assembly, is given by 



.(0) 



1 + 



1 + F{a,) (I 



„ 2d 



(33) 



We then compute the ZFC magnetization of an assem- 
bly of magnetic moments m,; = rii^B^i, i = 1,...,A/'. 
For this purpose we use the model introduced by Gittle- 
man et al. [76| which we summarize now. 

The ac susceptibility of an assembly of non-interacting 
particles, with a volume distribution and randomly dis- 
tributed easy axes, can be written as 



00 

x(T,a>)= j dVfiV) xv{T,L 



(34) 



where xo is the susceptibility at thermodynamic equilib- 
rium and xi is the initial susceptibility of particles in the 
blocked state [see [l| and many references therein]. The 
Fourier transform of psp leads to the complex suscepti- 
bility 



with the real part [76 



X 



(Xo + zcjrxi) 
1 -f iujT ' 



1 + uo^t"^ 



(36) 



(37) 



where ld is the angular frequency (= 2'kv). Starting 
from pT]) the application of an alternating field yields: 
a-) x' = Xo if ^ 1: At high temperature the magnetic 
moments orientate themselves on a great number of oc- 
casions during the time of a measurement, and thus the 
susceptibility is the superparamagnetic susceptibility xo- 
b) x' = Xi if 3> 1: At low temperature the energy 
supplied by the field is insufficient to reverse the mag- 
netic moments the time of a measurement. Here the sus- 
ceptibility is the static susceptibility xi- Between these 
two extremes there exists a maximum at the temperature 
Tlnax- x' can be calculated from (j37p using the formula 
for the relaxation time r appropriate to the anisotropy 
symmetry, and considering a particular volume V, one 
can determine the temperature Tmax- In an assembly 
of particles with a volume distribution, x' can be cal- 
culated by postulating that at a given temperature and 
given measuring time, certain particles are in the super- 
paramagnetic state and that the others are in the blocked 
state. The susceptibility is then given by the sum of the 
two contributions 



00 Va 

^[T,y)= j VVxi{T,V,iy) + J VVxo 
v„ 



,{T,V,iy), (38) 



where 14 = Vc{T,H) is the "critical volume" defined as 
the volume that discriminates between superparamag- 
netic particles of volume V < Vc and blocked particles 
with V > Vc, and is experiment-dependent. Equation 
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([38)) is then rewritten for the ZFC magnetization as fol- 
lows 



Vc 

M,faiH,T,^P) = JvV M,piH,T,V,^P) 



oo 

+ j VV Mb{H,T,V,iP) (39) 



where Msp and M;, are the contributions to the magneti- 
zation from the superparamagnetic and blocked particles, 
respectively, ip is the angle between the applied field and 
the anisotropy easy axis of each particle. 

The critical volume Vc can also be defined as the vol- 
ume at which the relaxation rate of the system is equal 
to the measuring frequency, i.e., F = In Ref. [73 we 
compute the contribution of DDI to Mgp and Mb and 
use the relaxation rate in p3[) to compute the DDI con- 
tribution to Vc- Then, we compute Tmax as a function 
of the applied field for increasing intensity of DDI, i.e., 
increasing DDI parameter S^d, within the limit of validity 
discussed above. 

In Fig. [T2| we present the results for an assembly of 
maghemite particles with mean diameter Dm = 5 nm 
and a distribution width 5 = 1.1. The curves are for 
a center-to-center inter-particlc distance that is a mul- 
tiple of D,n- The curve in closed circles is for the non- 
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FIG. 13: Magnetization of a polydisperse assembly of MSP 
particles with uniaxial anisotropy in the core (fee = Kc/J ~ 
0.0024) and TSA {k, = Ks/j '= 0.4). h = iiBH/J,t = 
kBT/J. The curve in up triangles is for A/" = 2551 and t = 
0.05. 



damping, as is the case in Fig. [T^l (A = 0.1), the prob- 
ability of switching increases. Then, if a magnetic field 
is added with increasing intensity the energy barrier is 
lowered and the magnetic moments switch at lower tem- 
peratures. 
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FIG. 12: The temperature rmax(J^) as a function of the ap- 
plied field for maghemite. 

effect of DDI is to change T^^xiH) from a bell-like curve 
with a maximum to a monotonic decreasing function, and 
this compares well with the experimental results [see Fig. 
1 of Ref. [73. As was stressed in Ref. 168, the effect of DDI 
is not only to change the potential energyscape, as was 
argued in many previous publications [l|, [80, [Slj , but also 
to introduce a transverse field that induces saddle points 
in the potential [79| . This in turn makes the relaxation 
rate quite sensitive to the damping strength, and for low 



In this last section we briefly comment on some recent 
(raw) results of the Monte Carlo calculations of the mag- 
netization of a polydisperse assembly of MSP nanopar- 
ticles as dealt with in section [IT| So far we have only 
considered non-interacting assemblies, and the work on 
the effect of DDI is still in progress. In Fig.[T3|we plot the 
(reduced) magnetization of two assemblies of 10 particles 
each. The total number of atomic spins in these assem- 
blies is 2551 and 4581. As the number of particles is the 
same in both assemblies, the particles in the assembly of 
2551 spins are smaller. These preliminary results show 
that the magnetization of the assembly with smaller par- 
ticles requires higher fields to saturate even at very low 
temperature, as can be seen by comparing the two up- 
per curves for t = 0.01 which corresponds to T = 1 K. 
As the temperature increases the spin disorder caused by 
the boundary and surface effects becomes stronger lead- 
ing to a decrease of the magnetization [see the curve at 
t = 0.05.1 



IV. CONCLUSIONS 

Contrary to the introduction, the conclusion is kept 
rather short. We will limit ourselves to emphasizing the 
unresolved issues owing to the complexities inherent to 
the systems under study, to the lack of sufficient expcri- 
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mental information, and to the deficiencies of the actual 
models. 

Nonetheless, we may contort ourselves saying that so 
far within the framework of the present models we have 
achieved a fair understanding of the static properties of 
nanomagnets, namely the spatial distribution of the mag- 
netization and its switching mechanisms. In particular, it 
becomes clearer, to both experimentalists and theorists, 
that the many-spin aspect of the nanomagnets is a pre- 
requisite for a better understanding of their properties 
when their size is reduced towards the nanometer. We 
believe that the work done so far by many authors has 
set the pace to future faster developments towards mas- 
tering and controlling the novel extraordinary features 
of nanomagnets in view of efficient and practical techno- 
logical applications. From the viewpoint of fundamental 
physics the facts established so far constitute a stepstone 
towards ever richer areas with more challenges to come. 

However, a great deal of work and endeavor still lies 
ahead of us because our actual knowledge of the subject 
is still plagued with quite a few obscure areas owing to 
the lack of sufficient pertinent information. This of course 
prevents us today from building a complete picture of the 
puzzle and makes our approach rather approximate. In- 
formation is gained by experimental investigation along 
with theoretical predictions and attempts of interpreta- 
tion. At present, the theoretical models developed so far 
present a few deficiencies. For instance, it is commonly 
assumed that the crystal structure on the surface is the 
same as in the core with the same atomic lattice param- 
eters. This cannot be wholly true considering the possi- 
bility of surface reconstructions. Of course, some models. 



including those presented here, do include apices, edges 
and facets, and the possibility of taking the exchange cou- 
pling on the surface as different from that in the core, or 
that between the core and the surface. However, there is 
no guarantee that it will become experimentally possible 
in the near future to estimate the atomic positions and 
lattice parameters and may be theoretically possible to 
perform crystal field calculations, and eventually check 
these assumptions. The magneto-crystalline anisotropy 
constant and exchange coupling in the core are also taken 
as those in the bulk underlying material. However, we 
have shown that the core of a nanoparticle does not en- 
joy the properties of the underlying bulk material. For a 
given material such parameters could vary with the ra- 
dius of the particle. In addition, the intensity and nature 
of surface anisotropy constitute a real challenge. There 
are many estimates but no firm understanding is achieved 
yet as to how surface anisotropy stems from the atomic 
structure in nanoparticles of reasonable size. 

For these interrogations to receive answers experimen- 
talists will have to devise ever more sensitive equipments 
and measurement techniques in order to probe the in- 
trinsic properties of nanomagnets and have direct access 
to the related physical observables. This, of course, will 
require new strategies for further isolating the nanomag- 
nets and rid them off the influence of their hosting matri- 
ces and mutual interactions. On the other hand, under- 
standing the influence on the particles static and dynamic 
behavior of the surrounding matrix and the intcr-particle 
interactions, is of paramount importance to efficient prac- 
tical applications. 
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